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ABSTRACT 

Context. The cosmological models for the formation of the first stars and the large scale structures now raise the question of how 
many dust particles were released to the general diffuse gas and how these impact the star formation process. In this framework, we 
focus on the scale of galaxy clusters. 

Aims. The aim of this study is to quantify the infrared luminosity of clusters as a function of redshift and compare this with the 
X-ray luminosity. This can potentially constrain the origin of the infrared emission to be intracluster dust and/or dust heated by star 
formation in the cluster galaxies. 

Methods. We perform a statistical analysis of a large sample of galaxy clusters selected from existing databases and catalogues. We 
coadd the infrared IRAS and X-ray RASS images in the direction of the selected clusters within successive redshift intervals up to 
z = 1. 

Results. We find that the total infrared luminosity is very high and on average 20 times higher than the X-ray luminosity. If all the 
infrared luminosity is to be attributed to emission from diffuse intracluster dust, then the IR to X-ray ratio implies a dust-to-gas mass 
abundance of 5 x 10~ 4 . However, the infrared luminosity shows a strong enhancement for 0.1 < z < 1, which cannot be attributed to 
cluster selection effects. We show that this enhancement is compatible with a star formation rate (SFR) in the member galaxies that is 
typical of the central Mpc of the Coma cluster at z = and evolves with the redshift as (1 + z) 5 . 

Conclusions. It is likely that most of the infrared luminosity that we measure is generated by the ongoing star formation in the member 
galaxies. From theoretical predictions calibrated on extinction measurements (dust mass abundance equal to 10~ 5 ), we expect only a 
minor contribution, of a few percent, from intracluster dust. 

Key words. Galaxies: clusters, Infrared: galaxies, X-rays: galaxies: clusters, intergalactic medium, diffuse radiation, Galaxies: star- 
burst 



(iDe Jonglll977l) . Finally, it has been conjectured since the ad- 
vent of molecular astrophysics dGould & Salpeterll 19631) that the 
grains surfaces play an important role in interstellar chemistry, 
and particularly in the synthesis of molecular hydrogen (see 
ICazaux & Tielensl d2002l) for a recent estimation). 

The question of whether dust could play a role in optically 
thinner media and on cosmological scales such as the inter- 
galactic medium (hereafter IGM) or the intracluster gas (here- 
after ICM), if present, has been s crutinized on a theoretical 
point o f view by severa l autho rs: see lDwek. Rephaeli & Mather I 
(119901). iPopescu et al.l (|2000|) . and IMontier & Giardl d2004l) . 
iMontier & Giardl (120041) computed the balance of dust heating 
and cooling with respect to the dust abundance: cooling domi- 
nates at high temperatures in the hot gas of virialized structures 
(i.e. galaxy clusters), and heating dominates in low temperature 
plasma under high radiation fluxes such as in the proximity of 
quasars. The details, of course, depend on the local physical 
parameters such as the grain size and gas density. The limita- 
tions of and uncertainties in the modeling are due to incomplete 
knowledge how both the formation of the first supernovae and 
the production of dust, correlate spatially and precede in time 
the heating of gas accreted onto the large-scale structures and 
dust destruction. The processes of feeding and enriching both 
the ICM and IGM with dust is even more of a key issue when 
taking into account that dust in these media is affected strongly 

by thermal sputtering and has lifetimes ranging from 10 6 yr to 

Send offprint requests to: M. Giard, e-mail: martin. giard@cesr . fr 10 9 yr (depending on the grain size). Therefore, when attempt- 



1. Introduction 

Dust is only a minor constituent of the known cosmic matter. 
It represents only 1% of the total dif fuse (non-stell ar) bary- 
onic mass in our Milky Way (see e.g. ISpitzeJ dl978l) . page 7 
and 8), and possibly 0.002% in the diffuse intergalactic medium 
dAguirre et alj 1200 ll) . However, small solid bodies (i.e. dust 
grains) are the main source of opacity for the electromagnetic 
radiation emitted by stars because they allow some electrons to 
be far more mobile than in free atomic or molecular orbitals. 
The practical consequence of this opacity is that dust thermal re- 
emission of absorbed photons has become a very efficient and 
popular tracer of the sta r form ation efficiency (see the review by 
Lagac heT Puget & Dold d2005[) and references therein). 

Moreover, the photodetachment of electrons from dust grains 
can occur more easily than from either a free atom or molecule, 
because the binding energies are smaller. This implies that 
the gas in which dust grains are embedded can be heated 
by photoelectrons if the stellar radiati on is sufficiently strong 
(see lWeingartner. Draine & Banj (12006) and references therein). 
Heating by the photoelectric effect of grains and cooling by 
their thermal infrared (IR hereafter) emission, are two opposite 
energy routes that, in addition to the gas processes, can have 
an important impact on the energy balance of gas clouds, and 
therefore on their ability to fragment, condense and form stars 
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ing to understand in detail the formation of large-scale struc- 
tures in the Universe, the impact of a dust component on the 
thermodynamics of the early ICM raises the important question 
of the heavy elements and dust grain synthesis in the first su- 
pernovae (SNe) and how they were able to seed the primordial 
diffuse gas and change its ability to fragment and form m ore 
stars and galaxies (see for instanc e ISchneider et al.l (120061) or 

iKampakoglou, Trotta & Siikl d2008l)). 

lElfgren & Desert! d2004l) and lElfgren. Desert & Guiderdonil 



(12001 " evaluated the dust emission in large scale structures from 
the time of reionization to the present by assuming the dust den- 
sity to be a fraction of the dark matter density and its mass dis- 
tribution to follow that of the dark matter, and by fixing the dust 
lifetime. Depending on the dust lifetime, these authors found that 
the integrated signal from the dust thermal emission produces 
brightness fluctuations on arcminute scales, which can be com- 
parable to CMB and galactic dust fluctuations. 

We therefore expect dust in the IGM/ICM to contribute to 
the diffuse IR emission in the direction of the large scale struc- 
tures, and particularly towards galaxy clusters. The only claim 
of detection of IR emission from intracluster dust was that 
of IStickel et all dl998l |2002|) towards the Coma cluster using 
the ISOPHOT instrument. A statistical detection of intracluster 
dust extinction towards S PSS quasars located behind clusters o f 
galaxies was obtained bv lChelouche. Koester & Bowed d2007l) . 
An upper limit was also derived bv lMuller et alj d2008l) by study- 
ing the statistical reddenin g of galaxies behind 458 RCS clusters. 
Finally. iBai et alj d2007bl) derived upper limits at 24 and 70 j-im 
from observations by the Spitzer telescope of the cluster Abell 
2029 (z = 0.08). Detecting emission from the ICM is indeed a 
difficult task because it must be distinguished from the IR emis- 
sion of the cluster galaxies. The overall cluster IR emission for a 
single cluster is comparable or lower than the average sky fluctu- 
ation produced by galactic cirrus clouds and the cosmic infrared 
background (due to background galaxies). 

From observational results, the IR emission from clusters is 
expected to originate pr imarily in star-forming cluster galaxies. 
iMontier & Giardl d2005l) performed a statistical detection of this 
overall dust emission in the direction of clusters by stacking the 
IRAS emission at 12, 25, 60, and 100 fim towards more than 
1 1 000 known galaxy clusters. They measured significant, but 
faint, emission at all IRAS wavelengths. This emission corre- 
sponded to the dust emission originating in both intracluster dust 
and dusty star-forming cluster galaxies. 

From measurements of the IR emission of star-forming 
galaxies f rom observa tions by the Spitzer telescope at 24 fim, 
IBai et al.l d2006l 120071) derived the luminosity function of in- 
frared galaxies in the Coma cluster (z = 0.02) and MS 1054-0321 
(z = 0.83). They demonstrated that the cluster galaxy luminos- 
ity function resembles that of field galaxies at both studied red- 
shifts. They also emphasized the strong increase in IR luminosity 
with redshift, which implied that the star formation rate (SFR) 
was higher in the past within cluster environments as it was for 
field galaxies. Their luminosity functions, when integrated over 
the entire range of galaxy masses, inferred a total infrared bolo- 
metric luminosity of 2.15 x 10 45 ft™ and 27 x 10 45 hzl erg /s re- 
spectively for Coma and MS 1054-0321 dBai et alj ( 120071) and 
private communication from the author). In contrast, the rest- 
frame X-ray bolometric luminosit ies of the two clusters ar e re- 
spectively 1.23 x 10 45 /z TO 2 erg/s dArnaud & Evrardlll999T) and 
3.05 x 10 45 h~l erg/s dKotov & VikhUninll2005h . 

In this paper we use the stacking method formerly developed 
bv lMontier & Giardl d2005l) to probe statistically the IR dust lu- 
minosity of galaxy clusters as a function of redshift. Our objec- 
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8 293 


4 725 
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Table 1. Selection of the cluster fields selected for stacking. 
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Fig. 1. Distribution in redshift of the stacked clusters for the 
three selected lists: NED (solid line), NSC (dashed-dotted libe), 
and SDSS (dotted line). See Sect. [2J] for the definition of the 
cluster lists. 



tive is to determine whether this luminosity evolves with red- 
shift and discriminate if possible between the contribution of the 
SFR in member galaxies and the diffuse emission from intraclus- 
ter dust. The stacked analysis was also completed on the X-ray 
Rosat All Sky Survey (RASS). We used the RASS as a reference, 
because the X-ray properties of galaxy clusters (i.e. structural 
and scali ng propertie s) ar e well studied a nd constrained (see re- 
views bv lVoitl d2005l) and lArnaudl d2005l) ). More specifically, it 
is well known that the X-ray lumin osity correlates w e ll wit h 
the temperature or mass of clusters dArnaud & Evrardl d 19991) . 
iMaughanl d2007l) . IVikhlinin et alj d2008l) or lPratt et alj d2008l) ). 

The paper is organized as follows. Sections [2] and [3]present 
the details of the data used, the selection of the fields stacked, 
and the stacking method. Our results are presented and discussed 
in Sects. |]and|5] 

In this paper, we use the following cosmology: Hq = 
70km/s/Mpc, Q,„ = 0.3 and Q A = 0.7. 



2. Data sets and catalogues 

2. 1 . The catalogues of clusters 

We used three different lists or catalogues of galaxy clusters as 
input to our stacking method. 

- The first one was a generic extraction from the NED 
databas43of objects classified as "clusters of galaxies". Each 
object has a coordinate and redshift. This list is extremely 



see the NED web site: http://nedwww.ipac.caltech.edu/index.html 
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inhomogenous. We started with a total of 18 829 references, 
which reduced to 13,746 clusters fields after source duplica- 
tion and overlapping fields had been identified (see Sect. [3]). 
Hereafter, we refer to this list as the NED sample. This is 
the same as the cluster l ist used in our previous investigation 
dMontier & Giardll2005l) . 

- The second selection corresponds to th e Northern Sky op- 
tical Cluster survey (NSC hereafter) by iGal et all d2003l) . It 
contains 8 155 references, each with an indicator of richness, 

N gals • 

- The third catalogue is the SDSS/MaxBCG cluster catalogue 
(SDSS hereafter) from lKoester et al.1 d2007l) . which provides 
a precise measurement of the cluster richness (N200, i.e. the 
number of galaxies within Rqw)- However, the redshift cov- 
erage of this selection is limited to [0.1, 0.3]. We begin with 
13 823 references. 

In Fig. Q] we provide the redshift histograms of three sets 
of stacked fields (see the description of the selection process in 
Sect.©. 

2.2. All sky survey data sets 

We performed our stacking analysis over five sets of all-sky 
survey data: the IRAS at 12, 25, 60, and 100 fim data and 
the ROSAT All Sky Survey broad band data (i.e. [0.1-2.4]keV 
RASS). For the IRAS survey, we used the reprocessed IRIS data 
set that was re calibrated with DIRBE and cleaned of the instru- 
mental stripes (Miville-Des cenes & Lagachell2005l) . 

For the RASS, we first derived count rate images that where 
computed to be the ratio of the count images (i.e. rs93*iml.fits) 
to the corresponding exposure map (i.e. rs93*mex. 

fitfl)- We then 

converted these count rate images into unabsorbed flux images 
by correcting for the absorption due to the hydrogen column den- 
sity in the [0.1-2.4]keV band. We used a fixed conversion factor 
from counts to fluxes of 1.5 x 10"erg crcT 2 for one count. This 
factor was deri ved using a WABS'-'MEKAL model with XSPEC (see 
lArnaudl 0996)) assuming a temperature of 1 .5 keV and a metal- 
licity of Z = 0.3 Z for the intracluster hot plasma at the median 
redshift of z = 0.2 (see Fig Q]), and a hydrogen column density 
of 3 x 10 20 cirT 2 , which was the average for the fields stacked. 

The value of the conversion factor was almost insensitive to 
the cluster redshift. From z = 0.01, there is a slight increase of 
~ 3% up to z = 0.4 and of ~ 10% up to z = 1. However, this 
is quite sensitive to both the intracluster gas temperature and the 
column density of the foreground HI. By varying the intraclus- 
ter gas temperature from 1 keV to 5 keV, the conversion factor 
is changed from 1.4 x 10~ n erg cirT 2 to 1.65 x 10~ n erg cm -2 , 
respectively. The most important effect originates in the fore- 
ground HI column density: the conversion factor rises from 
1.2xl0~ n erg cirT 2 to 1.8xl0~ n erg crrT 2 for Nhi that increases 
from 1.5 x 10 20 cnT 2 to 6 x 10 20 cirT 2 . 

Considering that the average HI column density of the fields 
stacked is 3xl0 20 cirT 2 with anrms fluctuation of 1.5xl0 20 cnT 2 , 
we became concerned about the errors introduced into the 
stacked fluxes when using a single conversion factor. We com- 
pleted the stacking analysis therefore using a conversion fac- 
tor for each cluster corresponding to its individual Nhi value. 
Surprisingly there was almost no detectable impact on the av- 
eraged fluxes obtained after stacking. The reason was that, al- 
though the conversion factor varies significantly from cluster to 

2 see the ROSAT Status Report #180 (ROSAT NEWS No. 70) 
on the web for a detailed description of the RASS products: 
http://heasarc.nasa.gov/listserv/rosnews/msg00126.html 



cluster, the averaging process is completed for such a high num- 
ber of clusters that the effect of any scatter on the averaged image 
diminishes to zero. The same was true for the cluster gas temper- 
ature, for which an initial scattering in the conversion factor was 
smaller. The average temperature used, 1.5 keV, was consistent 
with the average richness of the clusters that we selected from 
the SDSS entries, < /V200 > = 16, and the N20 to gas tem- 
perature correlation derived bv lRvkoff et al.ld2008l) for the same 
clusters (see their Table 4). 

After calculating the stacked fluxes, they were converted to 
luminosities in the [0.1-2.4]keV band using the associated lu- 
minosity distances computed for the adopted cosmology (see 
Sect. [Ql. 

Finally, we derived the X-ray bolometric luminosities from 
the broad band [0. l-2.4]keV luminosities using a correction fac- 
tor of 1.35, which was also computed by XSPEC using a MEKAL 
model for a cluster at z = 0.2 of temperature kT — 1.5 ke V and 
metallicity Z = 0.3 Z ). The bolometric correction was mostly 
sensitive to the cluster gas temperature with a conversion factor 
of 1.24 at 1 keV (or 2.07 at 5 keV). However, the significant size 
of our cluster sample enables us to use a single correction factor 
computed for the average gas temperature. 

We were able t o confirm our luminos ities by comparison 
with those given by iBohringer et al.1 d2004l) in the same energy 
band for clusters in the REFLEX catalogue; these numbered 
about one hundred of clusters with luminosities ranging from 
a few 10 42 erg/s to a few 10 45 erg/s. Our luminosities were sys- 
tematically slightly lower than the REFLEX ones with a mean 
difference of -16% and a rms scatter of +13%. 

3. Stacking the IR and X-ray maps in different 
redshit bins 

The entire stacking meth od used was similar to that described by 
iMontier & Giardl (120051) . We recall the main steps of the analy- 
sis: 

1. Distinct fields :We first eliminated multiple occurrences of 
the same field. This could occur for two reasons. First, be- 
cause several clusters were neighbors on the sky (at the cho- 
sen scale of our fields). Second, because the lists of clusters 
that we used, as in the case of the NED list, contained sev- 
eral references to the same object. In practice, once a field 
had been selected for stacking, it was identified by the coor- 
dinates of its center and any other field whose center coordi- 
nates were closer than 10 arcmin was not selected. 

2. Data missing : We then retained only fields for which data 
existed in all five sky surveys (i.e. the IRIS 12, 25, 60, and 
100 fim surveys and the RASS). This step eliminated 1 665 
fields from the NED list, 972 from the NSC catalogue, and 
1 545 from the SDSS catalogue. At this stage, we extracted 
a 1° x 1° sq. degree field for each of the relevant cluster di- 
rection. The pixel size of all extracted field was set to IRIS 
pixel size of 1.5 arc-minutes. The RASS images had an orig- 
inal pixel size of 0.75 arc-minute, which was rebinned to 1 .5 
arc-minutes by coaddition of pixels and nearest neighbor ap- 
proximation. This rebinning was performed to facilitate the 
data handling process for the different surveys. 

3. Galactic plane : We excluded all fields within a distance of 
10 degrees from the Galactic plane (NED references). 

4. Strong sources : In each field, we then masked all obvious 
sources in a 5 x 5 pixels square area. The sources were first 
taken from the IRAS point source catalogue. This was com- 
pleted by a source detection algorithm for both the IRIS and 
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Fig. 2. Radial profiles derived from the stacked images of the NED selection in three redshift bins centered at z = 0.09 (solid line), 
Z = 0.24 (dotted line) and z = 0.62 (dashed line) for the 100 fim IRAS survey (left panel) and the [0.1-2.4]keV RASS (right panel). 
The redshift bins were defined to each have equal numbers of galaxy clusters (about 830 - see text). The dashed vertical line shows 
the limit of the aperture used to compute the luminosities and the dotted-dashed horizontal line is the zero flux level. 



the X-ray images. This algorithm performed a 3x3 boxcar 
smoothing of each field and selected sources identified at a 
confidence level of 5 sigma with respect to the average. If 
an infrared source was detected at the center of a field (i.e. 
within a radius of 5 pixels) the field was removed from the 
list. However, we allowed fields to have an X-ray source de- 
tected at their center since massive galaxy clusters were well 
detected in the RASS data. This eliminated 1 651 more fields 
from the NED list, 947 from the NSC, and 1 063 from the 
SDSS lists. 

5. Background subtraction : At this stage, we performed a 
background subtraction for each map. We chose a 3 degree 
polynomial surface as background layer. The best polyno- 
mial was obtained by applying a least square fit to the image 
over an external ring which extended from pixel # 10 to pixel 
# 20 (i.e. between radii of 15' to 30' from center). The source 
mask defined previously in each sky direction was used to 
blank the image area used for the fit. The order chosen for 
the polynomial surface (3 rd order) is optimum because it is 
able to adapt to most backgrounds without altering data at 
the field center due to some uncontrolled oscillation in the 
fit. 

6. Noisy fields and weighting : Finally, we eliminated the nois- 
iest fields and weighted the remaining fields with their intrin- 
sic noise level at 100 fim wavelength. This was performed by 
computing the rms flux in each masked image over the outer 
ring defined for the background fit (between radii of 15' to 
30' from center). All fields for which this rms was higher 
than 5 times the standard deviation of all rms were then ex- 
cluded. This rms evaluation of the sky fluctuation performed 
on each field at 100 fim was also used to compute the weight 
of each image in the final stacking. For field#/ we used the 
following weighting: <x ; = {rms , + median{rmsi)y 2 . The me- 
dian of all rms values was added to each rms value to pre- 
vent any extremely flat images (which statistically happen) 
dominating the final sum. This final selection step eliminated 
2 120 more fields from the NED, 1 244 from the NSC, and 
1 898 from the SDSS lists. 

Our final data set contained 8 293 fields from the NED list, 
and respectively 4 725 and 7 476 from the NSC and the SDSS 



catalogues. The average Galactic HI column density in the di- 
rection of the selected NED clusters was about 2.8 x 10 2<) cirT 2 
(and 2.4 x 10 20 cnr 2 for the NSC and SDSS selection). 

Regarding the cluster richness, the average number of galax- 
ies in the NSC and SDSS clusters selected were respectively 
Ngais = 32 and A^oo = 16. From the correlation between A^oo 
an d the gas temperatur e derived for the SDSS cluster catalogue 
by iRvkoff et alj d2008l) (see their Table 4), we obtained a mean 
temperature of 1.5 keV for our SDSS selection. 

To study the evolution of the cluster luminosity with redshift, 
the selected cluster fields were sorted in order of increasing red- 
shifts. Redshift bins were then defined to have an equal number 
of clusters in each bin (that is 830 clusters in each of the 10 bins 
for the NED selection). The average flux image for each bin was 
calculated using the optimal weights specified above. The bins, 
weights, and source masks used in the stacking were the same at 
all wavelengths. 

When a positive signal was detected at the center of a stacked 
field at a given redshift, we computed the corresponding av- 
eraged luminosity. This luminosity was calculated by integrat- 
ing the flux inside the central 18' (9' radius from center) of 
the stacked images after subtraction of a possible residual back- 
ground computed to be the average flux in an outer annulus be- 
tween 9' and 18' radii from the field center. 

For the infrared data, the luminosity was computed for 
a wavelength interval of AA = A. This is a convenient 
quantity since astrophysical dust models show that the to- 
tal infrared luminosity can be well approximated as the sum 
of these contributions in the f our successive IRAS bands 
(iDesert, Boulanger & Pugetlll990l) . 



4. Results 

The stacking provided firm detections at all redshifts only in the 
60 fim, 100 fim and X-ray bands. 

Fig.|2]shows the radial profiles computed for the stacked im- 
ages of the NED selection at 100 fim and [0.1-2.4]keV for three 
redshift bins: z = 0.09, z = 0.24, and z = 0.62. There is a clear 
signal detected at all redshifts, at both infrared and X-ray ener- 
gies. 
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Fig. 3. Averaged stacked luminosities into redshift bins for the 
NED sample. Crosses are the IRIS 100 /mi AI^ luminosities 
and empty diamonds the IRIS 60 /urn ones. Stars correspond to 
the X-ray bolometric luminosities derived from the [0. 1 -2.4]keV 
RASS. The two curves show the total IR luminosity density 
(erg/s/Mpc 2 ) fro m the cluster galaxies c omputed with the SFR 
evolution law b y iLe Floc'h et alj d2005l) (dashed line) and by 
iBai et al.l d2007l) (dotted line). Both curves are normalized on 
the central luminosi ty density of the galaxy clusters Coma and 
MS 1054 as done in lBai et al.l d2007l) . 



In Fig. [3] we show the average infrared and X-ray luminosi- 
ties in each bin for the NED clusters sample. 

In Fig. 2] we show the ratio of the total infrared luminos- 
ity to the X-ray bolometric luminosi t y. Acc ording to the dust 
model of iDesert. Boulanger & Pugel dl990h . the total infrared 
luminosity (i.e. bolometric) was approximated by the sum of the 
60 and 100 /mi luminosity : Ljr =* Al, \(60 pmi) + AI,\( 100 pmi). 
This value should certainly be taken to be a lower limit for the 
total IR luminosity since we know from direct observation or 
from modelling that dust emits a large fraction of its luminos- 
ity longward of 100 ixm, both in norm al quiescent spiral galax- 
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the same plot indicate the IR-SFR to X-ray luminosity ratio de- 
rived by IBai et al.l d2007l) for two peculiar clusters, Coma and 
MS 1054. The infrared luminosity of these two clusters was de- 
rived from Spitzer 24 fim photometry of the member galaxies 
and does not include the diffuse emission if any is present. We 
assigned a 50% error bar to the IR luminosity of these two clus- 
ters as measured by the authors. 

4. 1 . Origin of the IR luminosity 

The first conclusion to draw from our results is that the average 
bolometric infrared luminosity of galaxy clusters is far higher 
than their X-ray luminosity. The average IR to X-ray ratio is 
about 20, possibly showing a slight increase with redshift (see 
Fig. Si. 

Our measurements of IR luminosity probe the contributions 
of all sources of infrared emission in clusters; i.e. star forma- 
tion in the member galaxies and any diffuse emission from intra- 
cluster dust. These two emission processes have different physi- 
cal origins and consequences. On one hand, infrared emission 
from galaxies is correlated with the total star formation rate, 
which provides a net energy contribution to the immediate en- 



Fig. 4. Ratio of the bolometric IR to X-ray luminosities in each 
redshift bin. The two stars indicate the I R-SFR to X-ray lumi- 
nosity ratio for two cluster as derived by IBai et al] d2007l) from 
Spitzer observations (Note that for these clusters the IR lumi- 
nosity is the integrated contribution from all infrared galaxies 
detected at 24fim). 



vironment. On the other hand, emission from intracluster dust is 
a cooling agent that can accelerate the condensation of diffuse 
baryons into stars and galaxies. 

High IR to X-ray luminosity rat ios were also m easured for 
two individual clusters of galaxies bv lBai et aLld2007l) : L;r/Lx — 
2 and 9, for Coma (z = 0.023) and MS 1054-0321 (z = 0.826) re- 
spectively. However, in these two cases, the IR luminosity is due 
only to the contributions of the individual cluster galaxies de- 
tected by the Spitzer-IRAC camera at 24|/m, and does not mea- 
sure any diffuse intracluster dust emission. 

Although this comparison relies on only two clusters, and 
should therefore be considered with care, it implies that the 
IR luminosity contribution of cluster member galaxies repre- 
sents a significant fraction of the IR luminosity measured in our 
stacked data, although possibly not all of the luminosity. Any 
additional IR emission could either be due to a population of 
unresolved galaxies in the Spitzer observations, s uch as IR emit- 
ting dwarf galaxies wi th low surface brightness dPopescu et al.l 
2002} Iffinz et al.ll2006h . or to diffuse intracluster dust. In the fol- 
lowing, we attempt to assess quantitatively the validity of these 
two hypotheses: the IR originates in (1) cluster galaxies, or (2) 
intracluster dust. 

We note that both Coma and MS 1054-0321 are well known 
non-relaxed clusters that show significant evidence of ongo ing 
dynamical activities dNeumann et a l. 2003; Gioia et aTll2004l) . 



4.2. Intracluster dust 

W e use the simple model of intracluster dust emission developed 
bv lMontier & Giardl d2004l) . where the dust grains are heated by 
collisions with free electrons and radiate in the infrared with 
properties typical of interstellar dust grains. To first order, the 
IR to X-ray luminosity ratio is proportional to Zj = Md ust jM gas , 
which is the dust-t o-gas mass ratio. I f we no rmalize on the curve 
shown in Fig. 8 of iMontier & Giardl d2004l) . the IR to X-ray lu- 
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minosity ratio over the relevant electron temperature range (10 7 
to IQ^keV) is then given by: 



tICM 

^ir 

tICM 



= 0.40 



(—) 

\ 10 5 / 



(1) 



To derive an IR to X-ray ratio similar to that measured, i.e. 
Ljk/Lx — 20, from the thermal emission of diffuse intraclus- 
ter dust alone, we therefore require a dust-to-gas mass ratio of 
5 x 10~ 4 . 

This dust-to-gas mass ratio would represent a significant 
fraction (~ 17%) of the total mass of metals available in an intr- 
acluster environment (assuming Z = 0.3 Z Q ). The precise dust- 
to-gas mass ratio in the intracluster medium cannot be easily 
constrained on purely theoretical grounds, since it is the result 
of a complex balance between the productio n and destruction 
proces ses of ICM dust (see the discussion in iMontier & Giardl 
(12004 V 

4.3. SFR of the member galaxies 

It is now well established that the infrared emission of galaxies 
increases significantly between redshift z = and z = 1, because 
the star formation rate (SFR) is measured to be higher in the 
past than the present. With Spitzer 24 pm data, the evolution of 
the infrared luminosity function of galaxies was characterized 
quantitat ively for both field gala xies in the Chandra Deep Field- 
So uth by |Le Floc'h et alj d2005l) and galaxies in massive clusters 
bv lBaiet all (120071) . 

We adopt the Schechter formalism dSchechterl 19761) to count 
the number of objects within each logarithmic (base ten) lumi- 
nosity interval: 

? ) (1 - ff) exp(-ZWL* ) 



<D /R (L) = <S>* IR {L 1R IL] R T- a > exp(-WL^) (2) 

where a — 1.41 is fixed as in iBai et alj d2006l) 

Bofh lLe Floc'h et all d2005l) and lBai et all d2007l) were able 
to describe the evolution of this luminosity function using power 



laws for L* lR and <b* lR : 



(l+zf (3) 

with a L = 3.2, a D = 0.7 inlLe Floc'h et alj d2005l) . or a L = 4.0, 
an — 1 .4 in lBai et all d2007l) . It is important to note that the two 
sets of coefficients are compatible within the reported error bars 
in these two works. 

We integrated the luminosity functions to derive the evolu- 
tion of the total SFR luminosity with redshift. The normalisation 
used at z — is adjusted from the so-called BvdHC sample of 
IR galaxies in t he Coma cluster (see filled circles in Fig. 5 of 
iBai et al.1 d2007l) ') : 

*(10 41 ers/ S> z = 0) = 59 /Mpc 2 / log 10 (L„) 

L* R (z = 0) = 10 10 - 49 L 

= 1.18 x 10 44 erg/s (4) 

The resulting t otal infrared SFR lum inosities are plotted in 
Fig, rj] as dashed dLe Floc'h et al.ll2005l) and dotted dBai et alj 
|2007|) lines. These integrated luminosities are actually close to 
a (1 + z) 5 evolution law. We emphasize that this normalisation 
of the SFR luminosity is relative to a projected cluster area of 
I Mpc 2 ; it is a coincidence that it corresponds to the same level 
as the total luminosities computed from our stacked IRAS im- 
ages at redshifts higher than 0.1. Although we do not exploit 
this coincidence, we believed that it was noteworthy for future 
quantitative investigations of this question. 



5. Discussion 

A clear interpretation of our results is difficult because the galaxy 
cluster sample that we extracted from the NED catalogue is obvi- 
ously affected by selection effects. For instance, the net decrease 
in the IR luminosities compared with the SFR model at redshifts 
lower than 0. 1 is likely to be a strong selection effect in the clus- 
ter mass function that produces a strong bias towards low mass 
clusters at low redshifts with respect to high mass clusters (i.e. 
the relative number of identified low mass with respect to high 
mass systems is far more important at low than at high redshift). 
This issue was the first motivation for stacking the cluster data 
at X-ray wavelengths and to use the X-ray data as a reference 
data set. In a second time, to tackle this problem further, we re- 
peated our stack analysis for the NSC catalogue, which has a 
more carefully defined selection function and provides an indi- 
cation of cluster richness: N ga i s . Finally, we checked with the 
SDSS cluster sample, which provides a more robust determina- 
tion of the cluster richness (i.e. N200), that N ga i s in the NSC is a 
valid measure of the cluster richness, at least for statistical stud- 
ies like this one. 



5.1. The L IR /L X and the dust-to-gas mass ratios 



As shown in Sect. I4.2I we require a dust-to-gas mass ratio of 
Zdust — 5 x 10~ 4 to explain the IR to X-ray luminosity ratio if 
dust in the ICM is the unique source of the IR emission. This 
ICM dust abundance is actually to be taken as an upper limit 
since we obviously know that there is also a contribution to the 
IR from the star formation in the member galaxies. The compar- 
ison with the IR-SFR to X-ray luminosit y ratios measure d for 
the two clusters Coma and MS 1054-0321 dBai et alj|2007h con- 
firms that the measured SFR may be insufficient to explain the 
high IR to X-ray luminosity ratio derived by our stacking anal- 
ysis. Moreover, our upper limit to the dust in the IC M is of the 
same order as that measured by iMuller et al.l d2008l) from their 
search for statistical reddening of background galaxies behind a 
sample of 458 RCD clusters with z < 0.5. Their extinction value 
converted into an upper limit to the dust mass within the central 
Mpc is 8 x 10 9 M . This value then translates into a dust mass 
abundance limit of Z,; < 2 x 10~ 4 . 

Conversely, one can use the existing constraints on the 
intracluster dust abundance to derive the Ljr/Lx ratio. By 
looking at the reddening of QSO s behind SDSS clusters, 
IChelouche. Koester & Bowenl d2007l) measured an average color 
excess of < E(g - z) >^ 3 x 10~ 3 mag. This detection is valid for 
the external parts of the SDSS clusters (i.e. R > /?200 - 1 Mpc). 
If we extrapolate this to the central Mpc of a cluster, this cor- 
responds to a total dust mass of M^ust = 3 x 10 8 M Q (see their 
Eq. 4). Compared with the corresponding gas mass in the same 
volume, this leads to Zj 10~ 5 . Note that our stacking inte- 
grates the luminosity within a disk area of 10' radius, which 
translates into 1 Mpc radius at z = 0.1 and respectively 4 Mpc 
radius at z = 0.6. The dust abundance derived from this extinc- 
tion measurement is consistent w ith the upper limit derived by 
iBovv. Hogg & Moustakas I J2008) from measurement of redden- 
ing of SDSS galaxie s behind galaxy cluste rs and with theoretical 
predictions given bv lPopescu et al.l (|2000|) . From Eq. (Q]), a dust 
abundance Z ( / 10~ 5 translates into j^- — 0.4. This is only about 
2% of the IR luminosity measured in the stacked data, leaving 
most of it (^ 98%) to SFR in the member galaxies. 
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Fig. 5. Same as Fig.[3]but for the NSC catalogue. 
5.2. Selection effects 

Obviously the cluster richness for the NED extraction increases 
with redshift as a result of ill-determined selection effects, and 
this will affect the stacked luminosities (IR and X-ray). 

To confirm this, we reproduced the entire stacking analysis 
for the NSC selection. The NSC spans a narrower redshift range 
than the NED selection (z < 0.4 for the NSC sample, and z < 1 
for the NED sample), but it includes an estimation of the cluster 
richness. To quantify the selection bias, we averaged the cluster 
richness in the same redshift bins as those fixed for the IR and X- 
ray luminosities. The stacked luminosities for the NSC sample 
are plotted in Fig. [5] We have retained the same redshift binning 
as for the NED analysis, meaning that we now have unequal 
numbers of clusters in each bin. We don't show the final two bins 
which have too few clusters to provide a statistical significant 
result (these two bins contain respectively 61 and 14 clusters, 
whereas we have more than 350 clusters in all other bins). The 
average richness in each bin, < N ga i s >, is shown in Fig. [6] 

The IR luminosities derived from the NSC selection show 
the same behaviour as those obtained from the NED selec- 
tion. The average richness in the bins indicates that the selec- 
tion is homogeneous for z > 0.1 with a mean NSC richness of 
< Ng a i s >- 35. As expected, there is a clear bias toward smaller 
clusters at z < 0.1, with a mean richness of < N g „i s >- 16 for 
the first bin, corresponding to z < 0.07. 

Finally, we checked that the richness parameter N ga i s given 
in the NSC catalogue is a robust probe of the cluster richness. 
We show in the lower panel of Fig. [7] that the X-ray luminos- 
ity correlates well with N ga i s . To compile this plot, we simply 
run the stacking again, but for clusters stacked into bins of equal 
N ga i s instead of redshift. However, for this analysis, the stacking 
procedure was modified to compute the background subtraction 
and luminosity of each cluster field individually, before stack- 
ing. The correlation between N ga i s and the infrared luminosity is 
good, but not as significant as for the X-ray luminosity, particu- 
larly because of the large error bar in the first richness bin. 

From this analysis, we derived the following power laws, 
correlating the X-ray and infrared luminosities to N ga i s : 

Lx = 10 42.6 ± 0.3 ( ^ ) 2.0 ± 0.4 £rg/s 

L IR = 1044.8*0.1(^0.7,0.2 erg/s (5) 



Fig. 6. Average richness, < N ga i s >, into the redshift bins for the 
NSC sample. 

We completed the same analysis for the SDSS selection and 
the corresponding cluster richness A^oo- The results are plotted 
in the upper panel of Fig. [7j They indicate that both the X-ray 
and the infrared luminosities correlate well with A^oo- This is 
not sur prising because the r ichnesses given for the SDSS clus- 
ters bv lKoester et alj d2007l) were carefully determined in a ho- 
mogenous way. The correlation analysis provides the following 
power laws: 



Lx = 10 43 - 3±01 (^) 1 - 5±0 ' 2 ^/, 

L/s = 10 44. 8± 0.06 ( A^0 )0 . 8± 0.2 e ^ /5 (6) 

If we compare our Lx/N^o relation with that derived by 
lRvkoffetal.1 00081) (their Fig. 6), we see that our results are 
compatible for A^oo - 50, given that our luminosities are bolo- 
metric and refer to h = 0.7 instead of h = 1 (i.e. L\, i = 
1.35L[o.i-2.4fev]/z~ 2 )- However, their slope is steeper because we 
did not scale our aperture in terms of either redshift or cluster 
richness in the same way as these authors ( 750 kpc or 7?2oo)- 
Therefore we measure slightly higher luminosities for the low 
richness clusters (A^oo < 30). 

The tight correlation between X-ray luminosity and richness 
is clear both for the NSC and the SDSS clusters and confirms 
that N g „i s in the NSC catalogue is a robust measure of the cluster 
richness, at least for statistical purposes as in Fig.|6]to investigate 
whether the cluster richness is biased with redshift. 

A secondary result of this analysis is that the exponent of 
the L/r to Ngais correlation is lower than for the X-ray (0.8 in- 
stead of 1.5 for the SDSS clusters). This is definitely related 
to differences in the physical emission processes that occur in 
the infrared and the X-rays. However, we do not investigate this 
property in this paper. 

6. Conclusion 

We have coadded IRAS and RASS images in the direction of 
thousands of galaxy clusters selected in the NED database for 
their low noise properties (smooth bac kground) following the 
same method as lMontier & Giardl d2005l) . We have compared the 
IR and X-ray luminosities, and investigated their evolution with 
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Fig. 7. IR and X-ray luminosities as a function of the cluster 
richness for the SDSS sample (upper panel) and the NSC sample 
(lower panel). The dotted and dashed-dotted lines show the best 
fit power law functions describing the data: the exponent is 1.5 + 
0.2 (SDSS) for the X-ray luminosities and 0.83 + 0.2 for the 
infrared. 



redshift. We showed that the averaged infrared luminosities are 
on average 20 times higher than the X-ray luminosities, and that 
the infrared luminosity increases with redshift approximately as 
(1 + z) 5 for redshifts between 0.1 and 1. We applied the same 
metho d to clusters selected from the NSC catalogue bv lGal et"ai1 
(120031) . This is a well selected catalogue of galaxy clusters for 
which a qualitative control of the selection bias with the redshift 
was possible. This showed the same behavior as the previous 
catalogue, but for a restricted redshift range 0. 1 < z < 0.4, which 
confirmed our result for the evolution of the IR luminosity. 

If all the infrared luminosity that we measure is attributed 
to emission from diffuse intracluster dust, this implies a dust-to- 
gas mass abundance in the intracluster gas of ^ 5 x 10 -4 . 
However, given other constraints on the abundance of intraclus- 
ter dust (Zd =a 10~ 5 from the exti nction towards the outer parts 
of the SDS S clusters measured by[C helouch eT Koester & Bowenl 
(12007b and lBovv, Hogg & Moustakasl d2008l) ). and the fact that 
the IR luminosity mimics the evolution of the SFR with red- 
shift, we believe that it is more likely that the main contribution 
to the cluster IR luminosity is from the ongoing star formation 
in the member galaxies. From theoretical predictions calibrated 
with the afore mentioned extinction measurement, we expect 



that only a minor contribution (a few percent) could originate 
in intracluster dust. This conclusion is supported independently 
by the overall shape of the spect ral energy distribution o f the 
stacked data that we published in Montier & Giardl (|2005l) , and 
which showed a high 12 to 100 /urn color ratio (of 0.062 for flux 
density, or 0.52 for luminosity). This high color ratio implies that 
nanoparticles emitting in the mid IR are abundant relative to big 
grains within the emitting phase. These tiny particles would not 
survive in hot intracluster gas, and are instead typical of normal 
spiral and star-forming galaxies. 
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